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ABSTRACT 



We report the discovery of WASP-Sb, a transiting planet of 2.25 ± 0.08 Mjup on a strongly inclined eccentric 8.15-day orbit, moving 
in a retrograde direction to the rotation of its late-G host star Evidence is found that the star is in a multiple stellar system with two 
other companions. The dynamical complexity of the system indicates that it may have experienced secular interactions such as the 
Kozai mechanism or a formation that differs from the "classical" disc-migration theory. 
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1. Introduction 

Transiting planets provide a wealth of information on the struc- 
ture and formation of planets. The measurement of planet ra- 
dius combined with its mass has found a surprising diversity 
in the mean densities and in particular "inflated" hot Jupiters. 
Spectroscopic measurement of the Rossiter-McLaughlin effect 
on the radial velocity during transits indicates that some of these 
planets may not be aligned with the rotation axes of their stars 
(see references in |Winn ( 2010[ l). The diversity in the observed 
spin-orbit misalignments is somewhat similar to that seen ear- 
lier in period and eccentricity distribution of planets detected 
by radial velocity surveys (see refs in |Udry & Santos (2007J 
and references therein). The recent sharp rise in the detections 
of transiting planets is the outcome of successful ground-based 
wide transit searches surveys among which WASP ( |Pollacco~et 
|al.|2006) is the most prolific. 



These discoveries have stimulated theoretical investigations 
of alternative formation scenarios to the migration theory (Lin 
et al.|199"6 Wu & Murray|2003) . These alternative theories ac- 
count for the discoveries of eccentric hot Jupiters on orbits not 
aligned with the rotation equator of their star (|Wu & Murray 
2003||Fabiycky & Winn 2009, Nagasawa et al.|2008[ [Barker & 



Ogilvie|2009| r 



Send offprint requests to: Didier.Queloz@unige.ch 

* Based on observations made with HARPSspectrograph on the 
3.6-m ESO telescope and the EULER Swiss telescope at La Silla 
Observatory, Chile 

** Radial velocity data are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-H A/| 



2. Observations 

2.1. The WASPS multiple stellar system 

The star WASPS (TYC2 7522-505-1) at q'(2000): 
23h 59m 36.07s, 5(2000): -35° 1' 52.9", was observed in 



2006 and 2007 by the WASP-south telescope ( |Pollacco et al. 
.2006| ). It is a y = 9.79 magnitude star with a Tycho (B-V) 
color of 0.73 which is indicative of a G8 spectral type. The 
Infra-red Flux Method (IRFM) (Blackwell & Shallis'l l977] l, 
using GALEX, TYCHO-2, USNO-Bl.O R-magnitude, and 
2MASS broad-band photometry, yields a distance of 87+7 pc. 

WASPS is identified in the CCDM catalogue (CCDM 
23596-3502) as the A component of a system of three stars. 
The B component is a 15th magnitude red star, 4 arcsec south 
of A, and the third component C is a 10th magnitude star (HIP 
1 18299, HD224664) 142 arcsec north of A. The radial velocity 
of HD224664 is 4.7 kms"' and stable over two years (Mayor 
priv com.) but differs from the WASP-8 value of -1.5 kms"'. 
The proper motions of the components also differs. It is there- 
fore unlikely that C and A are physically associated. 

We measured the photometry and position of WASPS and 
its nearby star (B component) with the Euler CCD camera of 
the 1.2m swiss Euler telescope at La Silla (see Fig. T]|. By com- 
paring with nearby stars, we obtained a magnitude difference 
Amy = 4.7, Am/ - 3.5. A separation and a projected angle 
was measured on the deconvolved images (Gillon et al. 2007[ l 
and we obtained 4.83 ± 0.01" and PA = 170.7 + 0.1°(only 
internal eiTors being considered). Assuming that WASPS and 
its B component are part of a multiple system, the color in- 
dices would represent those of an M star A similar photometric 
analysis of the individual 2MASS archive images indicates that 
Amy - 2.7, Aran = 2.2, and Am^: = 2.1, which are also indica- 
tive of an M star. The value mentioned in the Washington Visual 
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Fig. 1. Original (left) and deconvolved (right) V-band image 
from the Euler telescope of the A and B component of WASP-8. 



Double Star Catalog measured 70 years ago indicates 4.0" and 
PA=170°( |Mason et aLn2001) . This suggests little, if any, rela- 
tive motion of the two stars over the 70-year time span between 
these observations. When compared with the proper motion in 
right ascension of WASP-8, about 100 mas/yr ( [Zacharias et al.| 
2004 1, this indicates a common proper motion pair. Given the 
distance of WASP-S the sky-projected separation of the pair is 
about 440 AU. Using available differential photometry, we esti- 
mate the temperature of the B component to be about 3700 K. 



2.2. Photometric and radial-velocity observations 

WASPS was recorded simultaneously by two cameras of 
the WASP-south telescope during two seasons (2006 and 
2007). Altogether 11 224 independent photometric points were 
recorded with a typical sampling of 8 minutes. Transit events 
were detected in data from the first observation season. This 
triggered radial velocity follow-up observations of WASP-8 in 
November 2007 with the Coralie spectrograph mounted at the 
Swiss Euler tele scope ( |Baranne et al. 11996) |Queloz et al.|2000[ 



Pepe et al. 



2002| i. With a combined analysis of the radial- velocity 



data and the photometry including additional WASP data from 
the 2007 season, a transit period of 8.15 days was found. No 
changes to the spectroscopic profile were detected, ruling out 
a blended eclipsing binary or starspots as the cause of the radial- 
velocity variation ( [Queloz et al.||200T) (see bottom diagram in 
Fig |2]). In the next season, observations with CORALIE were 
continued, revealing an additional drift in the y velocity of the 
system (Fig.[3]l. No second-order curvature term was detected. 

On 25 August 2008, following up on the confirmation of 
the planet, a complete and densely sampled transit event was 
recorded in R band with the Euler telescope to improve the deter- 
mination of the transit parameters (Fig.|4i. On 4 October 2008, 
a spectroscopic transit was measured with the HARPS spectro- 
graph installed on the 3.6m telescope at La Silla. During the se- 
quence, 75 spectra (44 in the transit) were measured with an 
exposure time of 300 s, corresponding to a typical signal-to- 
noise ratio per pixel of 50. The radial velocity measurement 
from these spectra shows an obvious Rossiter-McLaughlin ef- 
fect with a shape suggesting a non-coplanar orbit (Fig. |4|. In 
addition, four spectra were measured on the same night before 
and after the transit to help us to determine the rate of change 
in the radial velocity outside the transit. Three measurements 
were obtained later at other phases of the system to improve 
the matching and zero point correction between CORALIE and 
HARPS data. During the measurement of the transit sequence, a 
significant change in telescope focus happened at JD 54744.592, 
improving the flux entering the fiber by a factor of 2. 




c»-0.04 
■^-0.06 

8.62 
5 8.60 
X 8.58 
g 8.56 
Li- 8.54 

8.52 



Fig. 2. Top: Overall Keplerian fit to the RV data for Wasp-8b. 
Black triangles indicates Coralie data and red dots the Harps 
data. The long-term drift was removed to plot the velocity 
in phase (the zero is set at the time of the transit). Bottom: 
Bissector span and FWHM (in km s 'unit) plotted with the 
phase of the orbit. The HARPS radial velocity data were shifted 
to correct from the y velocity difference with CORALIE . 
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Fig. 3. CORALIE radial velocity measurements (red dots) of 
WASP-8 superimposed on the best-fit solution (solid line). 



3. Determination of system parameters 

3.1. Spectral analysis 

The individual HARPS spectra were coadded and used for a de- 
tailed spectroscopic analysis of WASP-8. The results are dis- 
played in Table [T] As in previous WASP-papers ( Cameron et| 
al._2007) , the analysis was performed using the uclsyn spectral 
synthesis package and atlas9 models without convective over- 
shooting with Ha and //^ Na i D and Mg i b lines as diagnostics 
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Fig. 4. Top: Radial velocity measurement phased with the tran- 
sit (mid-transit is at 0). Black triangles are CORALIE data and 
red dots HARPS data. Bottom: Normalized transit photometry 
measurement of WASP-8. Black triangles indicates SuperWASP 
data and red dots the R-band Euler photometry data. The best-fit 
model is superimposed in blue. 



Table 1. Stellar parameters of WASPS derived from spectro- 
scopic analysis. The quoted error estimates include those given 
by the uncertainties in T^g, log g, and ^t, as well as atomic data 
uncertainties. 



Teff 


5600 ± 80 K 


[Na/H] 


-1-0.22 


±0.07 


logs 


4.5 ±0.1 


[Mg/H] 


+0.21 


±0.04 




1.1 ±0.1 kms-' 


[Si/H] 


+0.29 


±0.09 


V sin i 


2.0 ±0.6 kms"' 


[Ca/H] 


+0.24 


±0.12 


[Fe/H] 


-1-0.17 ± 0.07 


[Sc/H] 


+0.23 


±0.05 


log A(Li/H)+12 


1.5 ±0.1 


[Ti/H] 


+0.24 


±0.08 






[V/H] 


+0.30 


±0.08 


dist 


87±7 pc 


[Cr/H] 


+0.17 


±0.09 


age 


3-5 Gyr 


[Co/H] 


+0.29 


±0.07 






[Ni/H] 


+0.23 


±0.07 



1.2 
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Fig. 5. Comparison of the best-fitting stellar parameters from the 
transit profile and spectroscopic analysis with evolutionary mod- 
els interpolated at [Fe/H]=0.17. The isochrones are 100 Myr, 1 
Gyr, 5 Gyr and 10 Gyr. The evolutionary tracks are indicated for 
0.9, 1.0 and 1.1 Msun. 



for Teff and (log g). The abundances and the microturbulence 



were determined in a similarly way to the work of Gillon et al. 
(|2009| l and used as additional Tgff and log g diagnostics ( jSmalley 
|2005] i! 

The Li i 6708A line is detected in the spectra indicating an 
abundance of log A(Li/H) + 12 = 1.5 + 0. 1, which implies an 
age of 3-5Gyr accordi ng to |Sestito & Randich| ( |2005| l. However, 
IsraeUan et al. ( 2009[ ) noted that stars with planets have lower 



lithium abundances than normal solar-type stars, so the lithium 
abundance may not be a good age indicator for them. 

The rotational broadening v sin / was measured by fitting the 
observed HARPS profiles of several unblended Fe i lines. A typ- 
ical value of macroturbulence Vn,ac= 2 kms"' was adopted and 
an instrumental profile determined from telluric absorption lines. 
We found that v sin ; - 2.0 + 0.6 km s"' , which is typical of a G 
dwarf of intermediate age. 

3.2. Analysis of the planetary system 

This whole data set was found to detect without doubt a planet 
transiting the star WASP-8. We analyzed together the photomet- 
ric (WASP and Euler data) and the radial velocity data, includ- 
ing the spectroscopic transit sequence in this context. Our model 
was based on the transit modeling by |Mandel & Agof|p002| l and 



the radial velocity description by Gimenez (2006i. The best-fit 
model parameters and their error bars were computed using a 
MCMC convergence scheme that solves all parameters together. 
For details of the code and fitting techniques, we refer to Triaud 
let al.| ( |2009l l; |Collier Cameron et al.|p307 ) . To obtain a coherent 
solution, we determined the mass of the star by comparing the 
spectroscopically-determined effective temperature and the stel- 
lar density outcome of the MCMC adjustment, with evolution- 
ary tracks and isochrones of the observed metallicity from the 
stellar evolution model of Girardi et al. ( 2000 ). We converged 
iteratively on a stellar mass of 1 .04(+0.02 -0.09) M© and an age 
younger than 6 Gyr (see on Fig. |5]l. 

The free parameters of our model were the depth of transit 
D, the width of transit W, the impact parameter b, the period P, 
the epoch of transit centre To, the RV semi-amplitude K, e cos u) 
& e sin (D {e being the eccentricity and co the angle of the peri- 
astron), and V sin I cos /3 & V sin I sin /3, with V sin I being the 
projection of the stellar equatorial rotation, and /3 the projection 
of the angle between the stellar spin axis and the planetary or- 
bit axis. In addition, we employed free normalization factors for 
each lightcurve (WASP and Euler) and each set of radial velocity 
iyn for HARPS and yc for CORALIE ), which enabled variations 
to be made in instrumental zero points. From these parameters, 
physical parameters were derived to characterise the planetary 
system. The best-fit set of parameters that minimize the x^- (re- 
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Table 2. Fitted and physical parameters of the WASPS planetary 
system. The error bars are calculated at 68 % of the statistical 
distributions. 



D 


0.01276 


(+0.00033 -0.00030) 


W (days) 


0.1832 


( +0.0030 - 0.0024) 


b{R.) 


0.604 


(+0.043 -0.040) 


P (days) 


8.158715 


( +0.000016 - 0.000015) 


To (BJD-2450000) 


4679.33394 


( +0.00050 - 0.00043) 


/f (ms-') 


222.23 


( +0.84 - 0.60) 


dy/dt (ms"' yr"') 


58.1 


( +1.2 - 1.3) 


ecoso) 


0.02307 


(+0.0010 -0.0010) 


e sin CO 


-0.3092 


( +0.0024 - 0.0029) 


V sin I cos p 


-0.873 


( +0.059 - 0.064) 


rc" (ms-') 


-1565.76 


( +0.16 -0.21) 


7//''(ms"') 


-1548.10 


(+0.60 -0.13) 


V sin I sin/3 


1.59 


(+0.08 -0.09) 


RplRi. 


0.1130 


(+0.0015 -0.0013) 


RJa 


0.0549 


( +0.0024 - 0.0024) 


P* (Pe) 


1.22 


(+0.17 -0.15) 


R* (Re) 


0.945 


(+0.051 - 0.036) 


M, (Mo) 


1.030 


( +0.054 - 0.061) 


V sin 1 (km.s^^) 


1.59 


/ . f\ f\c\ f\ f\r\\ 

( +0.08 - 0.09) 


Rp/a 


0.00620 


( +0.00036 - 0.00033) 


Rp (Rj) 


1.038 


( +0.007 - 0.047) 


Mp (Mj) 


2.244 


( +0.079 - 0.093) 


a(AU) 


0.0801 


( +0.0014 - 0.0016) 


i(l 


88.55 


(+0.15 -0.17) 


e 


0.3100 


( +0.0029 - 0.0024) 


w(°) 


-85.73 


(+0.17 -0.18) 


)8(°) 


123.3 


( +4.4 - 3.4) 


computed at JDB= 


=2454691.15781 





duced is 0.86) are listed in Table [2] as well as their related 
computed physical parameters. With this best-fit solution one 
computes for the CORALIE data - 204 with 48 measure- 
ments, and for HARPS data - 188 with 82 measurements 
which implies that additional jittering is present that is not ac- 
counted for by the fitted model. Since the main deviation is re- 
lated to the CORALIE data, the uncertainties in the orbital so- 
lutions are most likely underestimated. However, the error bars 
in the Rossiter parameters are driven mostly by the HARPS on- 
transit data, and one can assume that they are almost correct. 

Our best-fit solution corresponds to a giant planet with an 
eccentric (e = 0.3) 8.16-day orbit and an additional long- 
term radial-velocity drift of 58 ms"' yr"'. The planet is dense 
with 2.25 Mj and a radius of 1. 04- Rj, in contrast to the sub- 
stantial fraction of 'inflated" hot Jupiters. Surprisingly, the pro- 
jected angle between the orbital and stellar spin axes is found to 
be y6 = 123.3°, indicative of a retrograde orbit. We note that 
VsinI = 1.59kms"' is in accordance with the line rotation 
broadening v sin / (in Table[T]) derived by the spectral analysis. 

We checked whether the partial defocusing of HARPS dur- 
ing the transit spectroscopic sequence had any effect on our re- 
sult. We divided the series into two subsets and considered for 
each of them an independent offset (y). We obtain a solution 
with a marginal improvement in the;\f^. By comparing the solu- 
tion obtained from these two sets with that for the complete set, 
the angle /3 was changed by 1 .5 cr. The defocusing problem does 
not affect the results of this paper. 



4. Discussion 

The detection of a hot Jupiter on an eccentric orbit that is mis- 
aligned with the stellar rotation axis and moving in a retrograde 



direction raises many questions about the origins of this sys- 
tem. Although the answer is beyond the scope of this paper, 
the visual faint companion and the drifting y velocity of the 
system are key components of the puzzle. From the observed 
separation between the A and B components, one can derive 
a most likely orbital semi-major axis (a = 1.35p ^ 600 AU) 
( [Duq uennoy & Ma yor|199l"| . The observed radial-velocity drift 
is therefore unlikely to be related to the B component of the bi- 
nary (y < GMa^^ < 1ms"' yr~'), suggesting that these is an 
additional closer companion of both unknown mass and period. 
The lack of curvature indicates that the companion is more mas- 
sive than the transiting planet. This intermediate body is very 
likely to play a significant dynamical role in the system. 

Apart from the complex dynamics of the whole system, the 
planet WASP-Sb is a "standard" hot Jupiter. It orbits a metal- 
rich star, which accounts for the observed increase in the inci- 
dences of hot Jupiters with the metallicity of the host star ( |Udry| 
& Santos|2007) . The period of WASP-8b is longer than the 3-4 
days typical value, but considering the eccentricity of its orbit, 
its periastron distance is typical of hot Jupiters. 

The orbit misalignment of the planet with the stellar rotation 
axis of WASPS is measured with the /? parameter The true angle 
between the axes of the stellar and planetary orbits is usually 
called i/' and is statistically related to p through sin / (unknown) 
and the orbital inclination (/) (see Fabrycky & Winn (2009') for 



details). When fi deviates significantly from zero, this provides 
us with a lower limit to the ijj. When (3 is beyond 90°, the orbital 
spin has the opposite direction to the stellar rotation provided 
that the orbit does not transit the star between its pole and its 
limb. According to Eq. 9 from Fabrycky & Winn ( 2009), this 
condition is met when / > 3.6 degree. By combining V sini, 
with the estimated age of the star, one can exclude such a small 
/ angle. Interpreted with the large /3 value, we can conclude that 
a true retrograde orbit is the most likely scenario for WASPSb. 

The origin of the unusual shape and orientation of the orbit 
of WASPSb is possibly related to the Kozai mechanism (Kozai] 
1962 Wu & Murray|20(j3) l or the outcome of a violent dynam- 
ical interaction history. The evidence of two other bodies and a 
possible series of secular effects (Takeda et al.1'2008 ) make the 
WASPS system unique and interesting for additional dynamical 
studies and a test case for formation scenarios of hot Jupiters that 
constitute an alternative to the disc-migration mechanism. 
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